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Abstract 

We describe impulsive gravitational pp-waves entirely in the distri- 
butional picture. Applying Colombeau's nonlinear framework of general- 
ized functions we handle the formally ill-defined products of distributions 
which enter the geodesic as well as the geodesic deviation equation. Us- 
ing a universal regularization procedure we explicitly derive regularization 
independent distributional limits. In the special case of impulsive plane 
waves we compare our results with the particle motion derived from the 
continuous form of the metric. 
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1 Introduction 

Plane fronted gravitational waves with parallel rays (pp-waves) are spacetimes 
admitting a covariantly constant null vector field, which can be used to write 
the metric tensor in the form jlj 

ds 2 = H(u,x,y)du 2 — dudv + dx 2 + dy 2 , (1) 

where u, v and x, y is a pair of null and transverse (Cartesian) coordinates 
respectively. In this work we shall deal especially with impulsive pp-waves which 
can be described by a profile function H proportional to a (^-distribution g], i.e. 
H(u, x, y) = f(x, y)S(u), where / is assumed to be smooth. Hence the spacetime 
is flat everywhere, exept for the null hypersurface u = 0, where it has a 8- 
shaped "shock". Such geometries arise physically as ultrarelativistic limits of 
boosted black hole spacetimes of the Kerr-Newman family ||, [Ij and multipole 
solutions of the Weyl family g . Moreover they play an important role in particle 
scattering at the Planck scale [|J. 

There are also intrinsic descriptions of impulsive pp-waves, most prominently, 
Penrose's "scissors and paste approach" ||, which consists in glueing together 
two copies of Minkowski space along the null hypersurface u = 0, identifying 
points according to (0, v, x, y) = (0, v+ (l/2)f(x,y),x,y). Penrose also intro- 
duced a different coordinate system in which the components of the metric 
tensor are actually continuous and, in the special case of a plane wave, i.e. 
f(x, y) = x 2 — y 2 take the form 

ds 2 = -dudV+(l-u+) 2 dX 2 + (1 + u+fdY 2 , (2) 

where u + — 6{u)u denotes the "kink-function" and is the step function. 
Clearly, the transformation relating (Q) with (||) has to be discontinuous, i.e. 
is given by 0| (for the case of a general pp-wave see 0]) 

x = (l + u+)X 

y = (l-u+)y (3) 
v = V + X 2 (u + l)8(u) + Y 2 (u- l)9(u) 

Hence -strictly speaking- the (topological) structure of the manifold is changed. 
In this work we begin our investigations using the original distributional form 
of the metric, motivated by the fact that physically, i.e. in the ultrarelativistic 
limit, the spacetime arises that way (cf. also the approaches of @, [|). Solving the 
geodesic as well as the geodesic deviation equation, we describe the geometry 
of impulsive pp-waves entirely in the distributional picture. As discussed in 
detail in jnj these equations, due to their nonlincarity and the presence of the 
Dirac <5-function in the spacetime metric, involve formally ill-defined products 
of distributions. To overcome this problem we use the setting of Colombeau 
algebras of generalized functions |0, [l2| which we briefly describe in Sec. ||. 
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Finally, for the special case of plane impulsive waves, we compare our results 
with calculations using the continuous coordinate system (see also [Q) to find 
that the particle motion in both spacetimes agree. Hence we argue that from a 
physical point of view the spacetime manifolds are identical. 

2 Colombeau Algebras 

In the vector space of distributions V no meaningful product can be defined. 
Moreover, as L. Schwartz showed in 1954, there does not even exist an 
associative and commutative differential algebra containing the space of ex- 
junctions as a subalgebra that allows a linear embedding of distributions. J.F. 
Colombeau |ll], [l2[ introduced differential algebras Q containing the space 
of distributions as a subspace, and the space of smooth functions as a faith- 
ful subalgebra. In the light of Schwartz's so called "impossibility result" this 
framework has to be regarded as the best one can hope for. 
To begin with, we give a short description of the algebra we are going to use in 
the sequel which provides a natural framework for studying nonlinear operations 
on singular data, hence singular (linear and nonlinear) partial differential equa- 
tions. The main idea is to regularize singular objects by sequences of smooth 
functions, i.e. in the space £ (fl) :— C oo (f2)' 0,1 ' ) , where denotes an arbitrary 
open set of R n . To achieve the consistency properties discussed above we intro- 
duce the spaces 

£ M (ty = {Me G (C oo (fi)) (0,1) : VK CCnVaeN^BN >0: 

sup | d a u e (x) \— 0(e~ N ) (e->0)}, (4) 

xEK 

Af(Q) = {(u e ) £ G (C* oo (^)) (0 ' 1) : VK CCQVaeN£,VM >0: 

sup I d a u e (x) |= 0{e M ) (e -► 0) } . (5) 

£jVf(f2) is a differential algebra with pointwise operations and Af(£l) is an ideal 
in it. We define the algebra of generalized functions, or Colombeau algebra, by 
the quotient 

0(fi) := £M(to)/tf(to) (6) 
and denote its elements by u = (u e ) e + Af(Cl) or-more carelessly-by u e . 
Distributions with compact support are now embedded into Q(R n ) by convolu- 
tion with a mollifier p e , defined as follows; let p £ £>(R") (Schwartz's space) 
with the properties / p(x) dx = 1 and J x a p(x) dx = Va G N™, |a|> 1, then 
we set p € {x) := (l/e n )p(x/e). So we have the map lq(ui) = (ui * p e ) t +7V(K"). 
This embedding can be "lifted" to an embedding i : D'(f2) Q(Q) by means 
of sheaf theory while smooth functions are embedded as constant sequences, i.e. 

<f) = (f)e- 

Next we briefly recall the concept of association in the Colombeau algebra, which 
is essential in applications. Two generalized functions u, v are called associated 
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with each other (u ks v), provided their difference converges to in T> . In 
particular, if u £ T>' then it is called the macroscopic aspect (or distributional 
shadow) of v. Equality in V is reflected as equality in the sense of association 
in Q, while equality in Q is a stricter concept (for example, all powers of the 
Heaviside function are distinct in the Colombeau algebra although they are 
associated with each other; moreover not every Colombeau function "casts a 
distributional shaddow," for example l(6) 2 ). 

The construction introduced above is known as the "special" or "simplified" 
variant of Colombeau's algebra. A more refined construction, i.e. the "full" al- 
gebra, avoids the dependence on p of the embedding of V . At the expense of 
some technicalities one can also construct a diffcomorphism invariant embed- 
ding. However, for the present application it will be sufficient to work in the 
"special" algebra. 

In the mathematical literatur Colombeau algebras have been extensively used 
to study nonlinear PDEs with singular data or coefficients (see eg. |l3| and 
the references therein). Recently there have also been a number of applications 
to general relativity, in particular the calculation of the curvature of cosmic 
strings p| |l7j , fL8| and the ultrarelativistic limits of Kerr-Newman black holes 
(see pj|, |20[ and references therein) . For a current overview of the topic see the 
review article of J. Vickers pl[ . 

3 Solving the Geodesic and Geodesic Deviation 
Equation in Q 

We start with a pp-metric of the form 

ds 2 = f( x l ) 6{u) du 2 ~dudv + {dx l f , (7) 

where x l [i — 1,2) denote the transverse coordinates. It is straightforward to 
derive the geodesic equation which (using u as an affine parameter, thereby 
excluding only trivial geodesies parallel to the shock) takes the form 

v(u) = f(x j (u))6(u) + 2dif{x j (u)) x l (u)S(u) , 
P(u) = ±dif(xt(u))8(u) , (8) 

where ' denotes the derivative with respect to u. 

Note that the first line involves a product of the <5-function with x l which we 
cannot even expect to be continuous (and in fact turns out to be proportional to 
the Heaviside function, cf. (l^]). We proceed by transferring the geodesic equa- 
tions into Colombeau's framework. The general strategy for solving differential 
equations in Q is to embed singularities (in our case: S) into Q, which amounts 
to a regularization, and then solve the corresponding regularized equations. In 
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order to obtain general results we are therefore interested in imposing as few 
restrictions as possible on the regularization of S. The largest "reasonable" class 
of smooth^] regularizations of <5 is given by nets (pe)ee(o,i) of smooth functions 
p £ satisfying: 

(a) supp(p £ ) {0} (e - 0) , 

(b) / p e (x) dx-^l (e ->• 0) and 




Pe{x) \ dx is bounded uniformly for small e 



(cf. the definition of strict delta nets in |p_3[] , ch. 2.7). Obviously any such net con- 
verges to 5 in distributions as e — > 0. To simplify notations it is often convenient 
to replace (a) by 

(a') supp(Pe) Q [-e,e] Vee(0,l). 

We shall call a net satisfying conditions (a'), (b) and (c) a generalized delta 
function. Denoting the (^-functions corresponding to x % and v by x\ and v t we 
are now prepared to state the following (cf. p2| ) 

Theorem 1 Let p f G Q(M) be a generalized delta function. The regularized 
geodesic equation 

v e (u) = f{x{(u)) p c {u) + 2dif(x{(u)) x\(u)p c (u) , 

xi(u) = ±d t f(xt(u))p e (u) (9) 

with initial conditions^ v e (— 1) = vq, v e (— 1) = i>o, 1) = x 1 q, 1) = £q 
/ias a unique locally bounded solution (v e ,X % e ) £ 5(K) 3 - Moreover these solutions 
satisfy the following association relations 

x\ w 4 + 4(i+w) + -9i/(4 + 4) u + 

v e « u + «o(l + w) + /(4 + 4)^(«)+ (io) 

+9i/(4 + 4) (4 + J^/(4 + 4)) u + • 

Hence from the distributional point of view the geodesies are given by refracted, 
broken straight lines. More precisely, the geodesies suffer a jump and a kink in 
the v- as well as a kink in the a^-direction when crossing the shock hypersurface. 

1 Note that, since T> is dense in L 1 , practically even discontinuous regularizations (eg. boxes) 
are included. 

2 Note that we have to impose initial conditions "long before" the shock. Geodesies starting 
at the shock cannot be treated in a regularization independent manner. 
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The scale of the effect is entirely determined by the values of the profile function 
and its first derivatives at the shock hypersurface and (in the special case of plane 
waves) reproduces exactly Penrose's junction conditions (|^). Note, however, that 
the above results are regularization independent even within the maximal class 
of regularizations of the Dirac 8. 

Our next goal is an analysis of the geodesic deviation equation for impulsive pp- 
waves in the framework of algebras of generalized functions. As in [[To) to keep 
formulas more transparent we make some simplifying assumptions concerning 
geometry (namely axisymmetry) and initial conditions. Writing x = x 1 and y = 
x 2 we suppose that / depends exclusively on the two-radius \J x 2 + y 2 and work 
within the hypersurface y = (corresponding to initial conditions j/o — — yo). 
Furthermore we demand vo — = xq. In this situation the Jacobi equation for 
the regularized deviation vector field N?(u) = (N?(u), N?(u), N* (u), N*(u)) 
takes the form 

K = N* = Nf = [N«f'(x e ) + \ N?f"(x e )} Pe + \ f'(x e )N?p e 

N: = 2lN?f(x e )p e y-N?f(x £ )p t + {N?f(x t )p e }~ (11) 
-N?f"(x e )x 2 p e - N?f'{x e )x ePe , 

where ' denotes derivatives with respect to r and we have suppressed the param- 
eter u. Furthermore x e is determined by (^|) with simplifiactions as discussed 
above. Note that the last equation in ( |TT| ) involves terms proprotional to 6 2 S 
and even S 2 . However, existence and uniqueness of solutions to the correspond- 
ing initial value problem in the Colombeau algebra as well as (regularization 
independent!) association relations (which we give explicitly for some special 
choices of initial conditions) are still guaranteed. 

Theorem 2 The Jacobi equation ( fl"^ j with initial conditions N^(—l) = n a and 
N?(-l) = n a has a unique solution in 0(R) 4 . If N?(~l) = and N?{-1) = 
(a, b, 0, 0) the unique solution satisfies the following association relations 

Nf « \af{x Q ){u + + e{u)) 

N: « b(l + u) + a{f(x Q )5(u) + ^f'(x ) 2 (8(u)+u + )}. (12) 

For initial conditions iV e a (— 1) = (0, 0, a, 0), N^(— 1) — we have 

Nf « a(l+«+) 

« af(x )(9+±f"(x )u + ). (13) 

Hence, viewed distributionally, in the first case the Jacobi field suffers a kink, 
a jump and a 5-like pulse in the u-direction as well as a kink and jump in 
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the x-direction overlapping the linear flat space behavior. These effects can 
be understood heuristically from the corresponding behavior of the geodesies, 
given by equation (fill). The constant factor a, which gives the "scale" of all 
the nonlinear effects, arises from the "time advance" of the "nearby" geodesies, 
represented by the initial velocity of the Jacobi field in the it-direction. However 
the second case shows that this is not the only effect producing kinks and jumps 
and we will refer to it in the discussion of the next Section. 



4 Discussion 

We have shown that the geometry of impulsive pp-waves can be described en- 
tirely within the distributional picture. Note that it was essential to use regular- 
ization techniques, i.e. Colombeau's algebra, to handle the ill-defined singular 
terms instead of introducing "ad-hoc" multiplication rules into Schwartz linear 
theory. Even within the maximal class of regularizations of the Dirac-5 we were 
able to derive regularization independent distributional geodesies and deviation 
fields. This (mathematical) feature also has a remarkable physical consequence. 
Interpreting the impulsive wave as a limiting case of sandwich waves of the form 
(cf. [p3[ ) H(u,x,y) = f(x,y)p e (u) — > f(x,y)S(u) we have shown that the im- 
pulsive limit is totally independent of the special form of the original profile p e 
(see also the results in Q). 

Finally (in the special case of plane waves) we discuss the relations of our ap- 
proach to the one using the continuous form of the metric. The metric (^]) has 
the advantage that simple particle motion can bee seen directly. Indeed free 
particles at fixed values of X, Y and Z = V — T after the shock start to move 
such that their relative X- and y-distance is given by the functions 1 + u + and 
1 — u + respectively. This is in total agreement with equation (|l3|). (Note that it 
is the coordinate transformation which introduces the motion in u-direction.) 
Moreover, we can solve the geodesic equations for the metric either by using 
the method of Sec. || or (since these equations only involve Heaviside and kink 
functions) by solving them separatly for u < and u > and joining them in 
a C 1 -manner. Either way leads to the solutions (u < 1, and using analogous 
initial values as before) 

X@ 

X(u) — x + xq(2 + u-) - 



1 + u, 



Y(u) = y + y U- + T ^— (14) 

V(u) = Vo + M l + u) + ^-^, 

1 — u 1+u 

where u_ := 6(—u)u. 

If we now formally transform equations (14) according to (|3|) we again obtain 



the geodesies (|11|) . Therefore we conclude that physically the two approaches to 
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impulsive plane waves, hence the two differential structures of the manifold, are 
equivalent. However, the transformation once more involves products of distri- 
butions ill-defined in the linear theory. Future work will be concerned with a 
mathematical analysis and interpretation of this discontinuous change of coor- 
dinates. Note that the transformation (|J) is given precisely by the limit of the 
geodesic equations (|ll| ) with vanishing initial velocities. Hence we can study the 
respective transformation in the sandwich case and then approach the impulsive 
limit. 
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